AD-A200  307 


BTiC  EILE  COEJ 


OFFICE  OF  NAVAL  RESEARCH 

Contract  N00014-83-K-0470-P00003 

R&T  Code  NR  33359-718 


Technical  Report  No.  129 

The  Behavior  of  Microdisks  and  Microring  Electrodes. 
Applications  of  Neumann’s  Integral  Theorem 

by 

M.  Fleischmann,  J.  Daschbach,  S.  Pons 


Prepared  for  publication  in  J.  Electroanal.  Chem 


Department  of  Chemistry 
University  of  Utah 
Salt  Lake  City,  UT  84112 


July  15,  1988 


DTIC 

SELECTE 
NOV  1  4 1988 


Reproduction  in  whole,  or  in  part,  is  permitted  for 
any  purpose  of  the  United  States  Government 


DISTRIBUTION  STATEMENT  A 

Approved  for  public  nltaM; 

Dlttrltmtion  Uulixnit<vf 


88  11  10 


iimMagBMaflaaaflagBH 


la  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  markings 


2b  DECLASSIFICATION/ DOWNGRADING  SCHEDULE 


NUMB£R<S) 


2  Distribution /availability  of  report 

Approved  for  public  release  and  sale. 
Distribution  unlimited. 


S  MONITORING  ORGANIZATION  REPORT  NUMBERISi 


6a  NAME  OF  PERFORMING  ORGANIZATION  6o  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 

University  of  Utah  (lt  *pplk*b>*) 


6c  aOORESS  (Cty.  Start,  amt  Z/P Cod*) 

Department  of  Chemistry 
Henry  Eyring  Building 


7b  AOORESS  (Cty,  Start,  and  /IP  Coot) 


8b  OFFICE  SYMBOL  I  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


(if  appAcaMt) 


kflElVlTltHnaiiKIfeH 


Ba.  NAME  OF  FUNOING/ SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


Be  AOORESS  (Cty,  Start,  and  ZIP  Code) 

Chemistry  Program,  Code  1113 
800  N.  Quincy  Street 


.TaKliMi*] 


’  Microring  Electrodes.  Applications  of  Neumann’s  Integral  Theorem 


»  M?W6»dfifl1k«^*.SDaschbach,  S.  Pons 


N000 14-83-K-0470-P00003 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT  I  TASK 

ELEMENT  NO  I  NO  I  NO 


WORK  UNIT 
ACCESSION  NO 


13a  type  OF  report 
Technical 


16.  SUPPLEMENTARY  NOTATION 


1 3b  TIME  COM89D  7/88  |14  DATE  dUl»ltS?Tl98gr,  Mont*.  Day)  hs.  PaCMcOuNT 

FROM _  TO _  I  I 


17 

COSATi  COOES 

FIELD 

GROUP 

SUBGROUP 

18  SUBJECT  TERMS  {Continue  on  reverte  if  neceua/y  and  identify  by  bloc* 

microelectrodes  ,  mass  transport 


19  ABSTRACT  [Continue  on  reverie  it  neceuery  and  identify  by  block  number) 

Attached . 


20  DISTRIBUTION/  AVAILAB>L:TY  of  ABSTRACT 
ED  UNCLASS. FiED/UNlIMITEO  □  SAME  AS  RPT 


22a  NAME  OF  RESPONSIBLE  NDiviDuau 

Stanley  Pons 


OD  FORM  1473,  3*  mar 


21  ABSTRACT  SECi.R.~y  ClaSSiF.CATiON 

□  ott  users  Unclassified 


SECURITY  CuASS  f  CA-  0* 


83  APR  ea  t  on  maj  oe  jiea  urr.  «»nausteo 
a, i  other  eait.ons  are  obsolete 


is  shown  that  Neumann’s  integral  theorem  can  be  used  to  derive  the 
steady-state  behavior  of  microdisk  electrodes.  Results  obtained  previously 
using  other  methods  for  uniform  surface  concentration  and  uniform  flux 
boundary  conditions  are  derived  as  limiting  cases  of  a  general  equation  which 
describes  irreversible  and  quasi-reversible  reaction  kinetics.  It  is  shown 


that  the  approach  can  be  used  to  assess  the  effects  of  the  tertiary  current 


distribution  and  to  predict  the  behavior  for  electrode  reactions  which  are 
non-linear  in  the  concentration  terms.  V'i.  O.  './yOT  w  '  I  *  -1  'Z  Z,  "tM  a  I  "C  '■ 

1  —  ■  1  1  V™  '  /  I  ^  1  1  \  i 

-  ■■■ 

\ 


Accession  For 

~stis  graai 

DTie  TAB 
Unanncuneed 
Ju3tlf lcatlcn— 


By - - - — 

Distribution/ _ _ _ 

Availability  Code9 
|Avali  and/or 
Dist  I  Special 


□  □ 


THE  BEHAVIOR  OF  MICRODISK  AND  MICRORING  ELECTRODES. 
APPLICATION  OF  NEUMANN'S  INTEGRAL  THEOREM  TO  THE 
PREDICTION  OF  THE  STEADY  STATE  RESPONSE  OF  MICRODISKS 


M.  Fleischmann 
Department  of  Chemistry 
The  University 
Southampton,  Hants.  S09  5NH 
ENGLAND 

John  Daschbacn,  and  Stanley  Pons 
Department  of  Chemistry 
University  of  Utah 
Salt  Lake  City,  UT  84112 
USA 


To  whom  correspondence  should  be  addressed. 


1 


Introduction 

— *  In  a  series  of  papers  we  have  recently  discussed  the  steady  state 

behavior  of  microdisk  and  microring  electrodes'^’  'as  well  as  the 

chronopotentiometric ,  chronoamperometric ,  and  A.C.  responses  of  microdisk 

electrodes'1  ’  and  the  determination  of  the  rates  of  homogeneous  reactions 

— ti 

(c.e.  mechanism)  coupled  to  electrode  processes  at  these  electrodes'  ’  .  Our 

approach  has  been  based  on  the  use  of  discontinuous  integrals  of  Bessel 

,  fTT"  {8  9) 

functions  (for  comparable  solutions  in  heat  conduction  see ~egv  ’  ')r  The 

steady  state  response,  for  example,  is  governed  by  the  solution  of  the 

diffusion  equation  in  cylindrical  coordinates 


D  +  5  ac  +  D  a!c  .  o 

dr2  r  dr  dz2 


[1] 


namely 


C  -  C*°-  J  f(a)  exp  (-az)  JQ(ar)  da  [2] 

J  o  y 

where  the  function  f(a)  is  chosen  to  satisfy  the  boundary  conditions.  The  two 
simplest  boundary  conditions  are  those  of  uniform  surface  concentration  over 
the  disk: 


C*  -  constant,  0<r<a,  z-0 


-  o 

dz 


[3] 

[4] 


r  >  a 


1 

I 

1 


4 


I 


4 


4 


2 


and  of  uniform  flux: 


D—  -  Q  (a  constant),  0<r<a,  z-0  [5] 

dz 


[2],  [3],  and  [4]  are  combined  with  the  discontinuous  integrals 


(or)  d a/a  - 


(aa>J  (or)  da 
1/2  o 


k/2 ,  r  <  a 

sin  1(a/r) ,  r  >  a 


[6] 

[7] 


associated  with 


r 


sin(aa) 


J  (or)  da 
o 

1/2 

J  (aa>J  (ar)  da 
1/2  0 


,  ..  2  2.1/2 

l/(a  -r  )  ,  r  <  a 

0,  r  >  a 


[8] 

[9] 


to  give 


*  2  .  ® 
C  -  C  -  -  (C 


f® 

-C*)  exp(-az)  sin(aa)J  (ar)d a/a 

*  n 


[10] 


Similarly,  [2],  [5]  and  [4]  are  combined  with  the  discontinuous  integral 


4 


4 


3 


*T> 

J  (ar)  J  (aa)  ua  - 
01 

0 


0, 

r>a 

[11] 

1 

[12] 

2a’ 

r-a 

1/a, 

r<a 

[13] 

to  give 


a> 

C  -  C*-  —  f  exp(-az)  J  (ar)  J. (aa)  ^ 
n  J  -  0  la 


[14] 


The  total  flux  to  the  surface  for  uniforn  concentration  conditions  can  be 
derived  from  [10]  as 


F 


4D(C®-CS)a 


[15] 


leading  to  the  definition  of  the  mass  transfer  coefficient 


(-.)«  -  iD/’*  I16> 

Similarly,  evaluation  of  the  average  surface  concentration  from  [10] 


C  -  c"-  8Qa/3wD 


[17] 


leads  to  the  definition  of  the  mass  transfer  coefficient  for  the  uniform  flux 


conditions 


r  x 


k* 


4 


-  3*D/8a 


[18] 


Mass  transfer  to  the  electrode  will  be  governed  by  [18]  for  irreversible  and 

quasi* reversible  reactions  at  low  overpotentials  and  by  [16]  for  all  reactions 

at  very  high  overpotentials.  As  we  have  indicated  elsewhere^ ,  the  effects 

of  the  distribution  of  potential  in  the  solution  (the  tertiary  current 

distribution)  will  cause  [18]  to  be  the  appropriate  choice  for  the  bulk  of  the 

polarization  curves  for  most  conditions;  the  effects  of  changing  the  nature  of 

the  assumptions  about  the  magnitude  of  the  mass  transfer  coefficient  are  small 

and  near  the  detection  limit  for  present  experiments^. 

The  approach  we  have  outlined  in  Equations  [2]  -  [18]  has  recently  been 
(23) 

extended  in  a  discussion  of  the  application  of  more  general  boundary 
conditions  such  as 

*The  discontinuous  integral  [11] -[13]  can  also  be  applied  to  the  solution  of 


ac  _  D  t  d  ac  t  p  a2c 

dt  dr2  r  3r  3z 2 


[19] 


for  a  range  of  non-steady  state  experiments  and  for  the  determination  of  the 

rates  of  homogeneous  reactions  in  solution  that  are  coupled  to  electrode 

reactions Chronoamperometric  and  chronopotentiometric  responses  for 

uniform  surface  concentration  conditions  have  been  investigated  using  a 

(10-22) 

variety  of  analytical  and  simulation  techniques  eg  .  The  results  of 

calculations  based  on  the  assumption  of  a  uniform  flux  over  the  surface  are 

again  close  to  those  bas«.d  on  the  uniform  surface  concentration  assumption  in 

(3  4) 

those  cases  where  comparisons  can  be  made  ’ 
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kC,  0<r<a,  z-0 


[20] 


to  the  interpretation  of  the  steady  state  behavior  of  disk  electrodes.  By 
analogy  to  [14]  it  was  assumed  in  that  work  that  the  weighting  function,  f(a), 
is  given  by  a  series  expansion 


J  (aa)  J  (aa)  J  (aa)  F'bJ  (aa) 

f(a)  -  bo  — |  —  +  b^—  +  b2-L—  +  \  -Till -  [21] 

(aa)  (aa)  (aa)  (aa)n  +  1 


(in  our  notation)  at  low  overpotentials  where  kinetic  control  will  be  dominant 
and  by  analogy  to  [10]  by  the  expansion 


f(a)  -  b 


J  (aa) 
1/2 

0  .  .1/2 
(aa) 


J  (aa) 
b  -3'2—  + 
1  (aa)3'2 


J»»(aa>  . ..  -V b-  J(H) 
2  <».)5'2  L 


(aa) 


[22] 


at  high  overpot er.tials  where  diffusion  will  be  dominant. 

It  is  evident  that  the  effect  of  boundary  conditions  more  complicated 
than  [3]  or  [5]  (eg  that  of  [20])  on  the  behavior  of  microelectrodes  (as  well 
as  the  influence  of  non-linearities  in  the  surface  reactions  including  that 
due  to  the  distribution  of  potential  in  the  solution)  can  be  readily  discussed 
in  terms  of  Neumann's  integral  theorem  (eg  see^  Chapter  14)  and  we  discuss 
the  application  of  this  theorem  to  the  prediction  of  the  steady-state  behavior 
of  microdisk  electrodes  in  this  paper. 


I h&.  Application  si  Neumann’s  Integral  Theorem 

The  behavior  of  microdisk  electrodes  in  the  steady- state  will  be  governed 


in  the  general  case  by  the  differential  equation  [1]  subject  to  the  boundary 
condition 


D~—  -  Q(p) ,  0  <  p  <  a,  2-0  [23] 

dz 

-  0.  p  >  a,  z  -  0  [2U] 

where  p  is  the  radial  position  in  the  plane  of  the  disk.  The  concentration 
generated  at  any  position  (r,z)  by  the  distribution  of  sources  Q(p)  can  be 
written  down  from  Neumann's  integral  as 


C(r,z)  -  [  ado  f* -SSP-L:.***!  3<£)  J  (ar)J  (ap)pdp  [25] 

*  o  *  o  ci  D 

where  integration  over  p  has  been  restricted  to  0  <  p  <  a  in  view  of  [24], 

The  derivation  of  the  complete  solution  therefore  depends  on  the  evaluation  of 

2 

Q(p)  using  the  appropriate  boundary  condition.  For  electrochemical 
experiments  we  are  concerned  with  a  distribution  of  sinks  over  the  surface  of 
the  disk;  a  convenient  distribution  of  such  sinks  can  be  based  on  the 
truncated  Fourier  series 


-Q(tf)  - 


a  cos(0) 


—  { 
sm  v 


c  +  c  cos(d)  +  c  cos(20)  +  •••  +  c  cos(n0), 
0  1  2  n 


[26] 


where 


p  -  a  sin(0) 


[27] 


Time -dependent  problems  can  be  discussed  in  terms  of  a  distribution  of 
sources  Q(p,s)  where  s  is  the  variable  of  the  Laplace  transformation  (for  an 


application  to  heat  conduction  see 


(24) 


> 


The  number  of  terms,  n+1,  in  the  series  is  chosen  so  that  C(r,z)  converges  to 
values  independent  of  n.  It  may  be  noted  that  a  choice  of  coefficients 


c  —  c  —  c  —  •••  —  c  ~  0 
12  3 


o 


[28] 


gives 


-Q(tf) 


-c 

0 


a  cos ( 0 ) 


2  2.1/2 
Ua  ~p  ) 


\ 


cf.  Equation  [ 8 ] J 


[29] 


[27]  and  [29]  in  [1]  give 


C(r.z) 


C 

J 

n  J  „ 


exp(-az)  JQ(ar)da 


-*l  2 

J  * 


a  sin(ff)  JQ(a  sin(0))dfl 


r 

*a 

0 

,2a, 

1  /  2 

exp(-az) 


J  (ar)J  (aa)da 
o  1/2 


[30] 


and,  with  [6]  and 


cQ  -  2D/w  (C®-C*) 


[31] 


we  recover  [10],  Similarly,  the  choice  of  the  coefficients 


I 


8 


c  -  c  -  c  -  -  c  -  0  [32] 

0  2  3  n 

gives 


-Q(*>  “  *c/a 


Thus  we  can  write 


1  ®  Ql 

C(r ,z)  -  -  -  J  exp(-az)  JQ(ar)da  J  c 


c  .<*> 

—  exp(-az)  J  (ar)J  (aa) 
D  Jo  oi 


da 

a 


and,  with  [13]  and 


[33] 


[34] 


cx  -  Qa  [35] 

( 

we  recover  [14],  Thus  the  first  two  terms  of  the  series  [26]  give 
respectively  the  diffusion  controlled  limit  at  high  overpotentials  and  the 
4  kinetically  controlled  limit  at  low  overpotentials.  It  would  be  expected, 

therefore,  that  a  relatively  short  truncated  Fourier  series  will  adequately 
represent  the  behavior  at  intermediate  overpotentials. 


At  these  intermediate  overpotentials  we  need  to  evaluate  the  integral 


9 


f  a  sin(0)J  (aa  sin(0))/c  +  c  cos (0)  +  c  cos(20)  +  •••  +  c  cos(n0)ld0 

J  o  l  0  1  2  n  J 


It  is  convenient  to  express  the  terms  c^costjfl)  in  this  series  in  terms  of 
powers  of  cos(0)  using 


cos(j0)  -  Real  part  of  (exp(ij0))  -  Real  part  of  (cos(0)  +  i  sin(0))J 


and  Q(q)  can  be  reexpressed  as  a  series  in  cos^q 


-Q(9)  - 


1-  { 

os  (6)  ' 


a  cos(0) 


d  +  d  cos(0)  +  d  cos2(0)  +  •••  +  d  cosn(0) 
0  1  '  2  n 


where  the  coefficients  d^  are  expressed  in  terms  of  c^ .  The  nature  of  the 

coefficients  d^  therefore  depends  on  the  length  of  the  Fourier  series  [26] 

chosen  to  represent  Q (0);  values  of  the  coefficients  for  n  up  to  10  are  listed 

in  Table  1.  The  advantages  of  using  the  particular  series  expansion  [38]  to 

/  26  \ 

represent  the  flux  will  be  discussed  elsewhere' 

The  Neumann  integral  can  be  rewritten  using  [38]  giving 


2 

C(r,2)  -  -  —  exp(-az)J  (ar)da 
D  "  o  0 


£i? 


Q(0)Jfl(aa  sin(0))sin(0)cos’i(0)cos(0)d0 


$ 


D  -1  o 


J  exp(-az)Jo(ar)da  J  J  d^J^ 


aa  sin(0))sin(0)cos  (0)d0  [39] 


Using  the  known  result  (compare  eg 


(7,25) 


-*/  2 

J  JQ(oa  sin(*))sin(*)cosJ(*)d*  - 


2<(j-x)/2>r((1+1)/2) 
/  * ( ( 3*1 > / 2  ) 


we  derive  the  values  of  the  integrals  listed  in  Table  2;  we  express  these 


terms  as 


h 


(J+1J/2 


and  note  that  these  values  are 


independent  of  the  length  of  the  Fourier  series.  [26].  We  combine  the 
appropriate  choice  of  the  coefficients  in  Table  1  with  the  coefficients  in 
Table  2  to  give 


C(r ,z) 


-it  l>. 


J r  \(aa) 

f— 

exp(-az)Jo(or)  1 
(aa) 


In  the  application  of  [41]  we  are  especially  interested  in  the  values  of  the 
concentration  at  the  surface  of  the  electrode  where 


The  Ueber-Schaftheitlin-Sonine  type  integrals  can  be  conveniently  evaluated 

(7  25)  (2) 

in  terms  of  hypergeometric  (Gauss)  series  ’  (compare  )  using 


J  f  \ (oa) 
a  (j*l> ' 

2  1/2 

J  (or) - - da  -  - - -  2 

0  (a.)<(JM,/2)  2,<J*l,/2,r((J+2)/2)a 


'&i'  -  3 


As  can  be  seen,  the  value  of  any  particular  integral  at  a  position  r/a  depends 
on  j  and  we  denote  this  value  by  f^(r/a),  see  Table  2.  Ve  can  therefore  write 
[42]  as 


D  l__t 


n  \ 

1/2  V'  T  —  .  2x 

-  -  - -  \  d  F  -i;  1;  -  1 

2D  J  r  2  H2  2  a2/ 


Equation  [44]  provides  the  link  between  the  assumed  form  of  the  flux  [26]  (or 
[38])  and  the  concentration.  In  the  case  of  relatively  simple  boundary 
conditions  such  as  [20],  evaluation  of  the  concentration  at  (n+1)  radial 
positions  denoted  by  rk  provides  the  (n+1)  linear  equations  required  to 
determine  the  (n+1)  coefficients  (d  ...d  )  or  (c  .  ..c  ).  For  simplicity,  we 


2 


write  the  set  of  equations  in  term  of  the  coefficients  d 


ac 


v3zy 


“  Q(*k)  - 

r-r 


a  cos(0k) 


d  +  d  cos(0  )  +  •••  +  d  cos  (8  ) 
0  1k  ok 


M 


1/2 


2  2 
a  -r 

1/2 

2  2 
a  *r 

A  4.  A 

k 

+  d 

2 

k 

U  T  U 

0  1 

2 

a 

2 

a 

. 

.  4 

+  •  •  •  +  d 


2  2 

A  -  Y* 

n/2- 

a  i 

k 

2 

a 

.  , 

-  kC 


[45] 


j-o 


The  coefficients  d^(or  c^)  can  therefore  be  determined  by  a  set  of  equations 
such  as  [45]  and  the  behavior  of  the  system  is  then  fully  specified.  Thus  the 
local  flux  is  given  by  any  of  the  equivalent  forms  [26],  [38],  or 


2  2^ 
a  -r 

kJ 


1/2 


2  2 
a  -r 

1/2 

2  2’ 
a  -r 

H  4-  H 

+  d 

2 

U  T  u 

0  1 

2 

a 

2 

a 

. 

.  . 

.  j 

+  •  •  •  +  d 


2  2 
a  -r 


n/2. 


[46] 


while  the  total  flux  is  given  by  any  of  the  relations 


13 


d  m 

' 3z ' 

n 


1C- 


2*D  —  rdr  -  2*Da  )  c  cos(j$)sin(0)d0 


V  f>*  > 

a  )  d  cos 

U>  J 


(«)sin(*)d* 


°  (  „  ,\1/2 

D>M  F, 


-  2ira 


■I* 


The  concentration  distribution  C(r,z)  is  given  by  an  equation  such  as  [41] 
while  the  surface  concentration  is  given  by 


C(r)  -  C 


D 


(or  by  equivalent  forms  written  in  terms  of  6). 


I 


More  Complex  Boundary  Conditions 
Quasi -reversible  Reactions 

Certain  examples  of  the  application  of  more  complex  boundary  conditions 
can  be  analyzed  by  closely  similar  methods  to  those  outlined  above.  For 
example,  for  a  redox  couple  following  quasi -reversible  Butler-Volmer  kinetics 


[49] 


we  obtain 


MU  -  “-(?)[■  -  “'01 


— —  +  -5— 

DC*  D  C\ 

0  0  R  R' 


[50] 


vj  -o 


We  therefore  again  obtain  a  set  of  (n+1)  linear  equations  which  can  be  solved 

for  the  coefficients  d  (or  c  ).  Once  these  coefficients  have  been 

i  J 

determined,  the  concentration  distributions  and  total  current  at  any  given 
overpotential  follow  as  above  from  [41],  [48],  and  [47]. 


Xh£  Tertiary  Current 


Analysis  of  the  behavior  of  microdisk  electrodes  based  on  the  solution  of 


the  diffusion  equation  above  (i.e.  neglecting  the  effects  of  the  distribution 
of  potential  in  the  solution)  lead  to  the  conclusion  that  the  flux  at  the 


15 


edge  of  Che  disk  becomes  infinite  in  Che  limiting  current  region,  cf .  Equation 

[8].  Such  infinities  are  impossible,  however,  as  the  local  overpotential 

would  then  also  have  to  become  infinite.  In  reality,  the  "throwing  power”  of 

the  system  in  this  potential  region  must  be  sufficiently  high  to  make  the 

current  distribution  much  more  uniform  than  that  given  by  [8].  We  have 

pointed  out  chat  a  consequence  of  the  tertiary  current  distribution  therefore 

is  that  the  solutions  based  on  a  model  of  a  uniform  flux  over  the  surface, 

equations  [14],  [17],  and  [18],  is  more  likely  to  apply  over  the  bulk  of  the 

polarization  curve  than  solutions  based  on  a  modal  of  a  uniform  surface 

(3) 

concentration  over  the  surface.  Equations  [10]  and  [16]  .  It  is  clear  that 

an  adequate  analysis  of  the  behavior  of  microdisk  electrodes  at  high 
overpotentials  must  take  due  account  of  the  tertiary  current  distribution. 

As  has  been  shown,  the  use  of  the  Neumann  Integral  [25]  coupled  to  a 
suitable  Fourier  series  expansion  of  the  local  flux  [26]  leads  to  a  formal 
solution  [41]  or  [42]  which  links  the  distribution  of  the  flux  over  the 
surface  to  the  local  concentration.  The  coefficients  in  the  Fourier  expansion 
can  therefore  be  determined  by  applying  the  relevant  boundary  conditions  at  a 
sufficient  number  of  positions  on  the  surface.  This  procedure  can  evidently 
be  extended  to  take  into  account  the  effects  of  the  distribution  of  potential 
in  the  solution  especially  that  at  the  surface  of  the  disk. 

In  the  presence  of  excess  support  electrolyte,  the  potential  in  the 
solution  will  be  governed  by 


dr2  r  dr  dz2 


[51] 


which  is  entirely  analogous  to  the  differential  equation  [1]  governing 
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r 


r. 


diffusion.  We  take 


<t>  -  0,  r  -  z  -  ® 


[52] 


and  note  that  the  local  current  density  at  the  surface  of  the  disk  can  be 
written 


I  -  zFQ(p) 


[53] 


The  potential  at  any  position  in  the  solution  can  therefore  be  written  as 


*(r,z)  *  Jada  Q<*)J#<ar)J#<ap)*d/> 


[54] 


o  o  a 


and  at  the  surface  of  the  disk  we  have 


O)  £ 

^(r)  -  ~  |  JQ (or)da  f  Q(p)Jo(ap) pdp 

"  *  a 


:s5j 


Again,  using  the  representation  of  the  local  flux  by  [26]  and  [28]  we  obtain 
by  the  same  methods  as  those  used  to  discuss  the  concentration  distribution 
[48] 


56 


*  £  i  J 

j-o 


In  this  paper  we  will  restrict  attention  to  the  generalization  of  the  boundary 
condition  [20]  to  take  account  of  the  potential  distribution  [56].  Ve  now 
regard  the  rate  constant  k  as  the  appropriate  value  uniform  over  the  surface 
of  the  disk  provided  we  can  neglect  the  potential  distribution  in  the 
solution,  l.e. 


k  -  k  exp f-Ptf | ,  ^(r)  -  0,  0  <  r  <  a 

'  RT  ' 

For  such  a  cathodic  Tafel-type  relation  we  can  therefore  write 


[57] 


k(r)  -  k  expf~a^-r)F°| 
l  RT  7 


-  k  exp 


2 

-azF  a 
*cRT 


n 


J-0 


[58] 


As  before,  the  (n+1)  coefficients  of  the  Fourier  series  must  be  determined  by 
the  application  of  the  boundary  condition  [20]  in  this  case  modified  by  [58]. 
For  example  we  can  use 


r 


V 


. 


L 
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'  ' 

dc 

1 

* 

2  2 
a  -r 

1/2 

2  2 
a  -r 

2  2 
a  -r 

n/2- 

(*Hf! 

d  +  d 

0  1 

k 

+  d 

2 

E 

k 

r-r 

k 

2 

a 

.  4 

2 

a 

k  i 

T  •  •  *  T  O 

n 

2 

a 

.  j 

-  W 


exp 


j-o 


-°zF  -  Y^d  e  f  f— )  .  .  .  c"  -  -  V"d  e  f  f— ) 
icRT  J  J  H*JJ  J  J  n«J 


J*o 


[59] 


For  small  perturbations  of  the  bulk  concentration  and  small  values  of  the 
potential  in  solution  (i.e.  for  the  lov  overpotential  region)  sets  of 
equations  such  as  [59]  can  be  linearized  and  solutions  for  the  coefficients  d^ 
can  therefore  again  be  obtained  by  matrix  methods.  At  high  overpotentials  it 
is  however  necessary  to  solve  the  set  of  non-linear  equations.  This  can  be 
done  iteratively,  the  set  of  the  coefficients  d^  in  the  argument  of  the 
exponential  function  being  determined  by  the  previous  cycle  of  the 
calculation;  the  same  matrix  methods  can  therefore  still  be  applied.  Values 
of  the  concentration  and  total  current  are  derived  as  in  the  other  cases  when 
convergence  of  the  coefficients  has  been  achieved;  equally,  potential 
distributions,  average  potential  drops,  etc.  can  be  predicted  using  [56]. 


1 


Boundary  Conditions  Non-linear  in  the  Concentration  Terms 

The  kinetics  of  electrode  reactions  are  frequently  non-linear  in  the 
concentration  terms  either  because  of  the  complexities  of  the  reaction 
mechanism  (total  order  of  reaction  *  1)  or  because  of  adsorption  of  the 
reactants.  The  strategy  outlined  in  the  previous  section  can  also  be  used  to 
derive  complete  solutions  in  these  cases  and  we  illustrate  the  approach  for 
the  simple  example 


J 
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-  kC2  -  kC'  •  C  [60] 

dz 


We  evaluate  the  set  of  equations  linear  In  d^ 


ac) 

1 

• 

2  2 
a  -r 

1/2 

2  2 
a  -r 

2  2  I 

a  -r 

n/l * 

1  ™ 

:)*'* 

d  +  d 

0  1 

k 

+  d 

2 

ft 

+• • •+  d 

s 

k 

dz 

i  r-r 

v 

2 

a 

2 

a 

2 

a 

ft 

» 

.  j 

.  . 

.  , 

. 

[61] 


where  the  concentration  C'  and  the  coefficients  d'  have  been  evaluated  in  the 


previous  cycle  of  the  calculation. 


J 


It  is  apparent  that  the  use  of  Neumann's  integral  allows  a  detailed 
analysis  of  the  behavior  of  microdisk  electrodes.  The  choice  of  the  Fourier 
series  [26]  (or  the  equivalent  forms  [38]  or  [45])  to  express  the  local 
reaction  rate  appears  to  offer  special  advantages:  the  results  for  the  uniform 
surface  concentration  [3]  and  uniform  surface  flux  [5]  boundary  conditions 
follow  as  special  cases  and  a  single  series  such  as  [41]  can  be  seen  to  give 
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Che  concentration  distribution  under  all  conditions  (rather  than  the  two 
series  [21]  and  [22]  which  have  been  applied  to  the  low  and  high  overpotential 
regions,  respectively^^ ) .  ^ 

A  particular  advantage  of  this  approach  is  that  it  allows  a  systematic 
analysis  of  boundary  conditions  non-linear  in  the  concentration  terms  as  well 
as  of  the  effects  of  the  tertiary  current  distribution.  This  latter  analysis 
gives  a  precise  description  of  the  ohmic  potential  drops  in  solution:  these 
ohmic  potential  losses  (eg.  that  to  the  edge  of  the  disk)  will  be  quite 
different  to  those  at  microspherical  electrodes  and  the  magnitudes  depend  on 
the  disk  radii.  Such  predictions  are  in  line  with  experimental  observations 
but  are  contrary  to  currently  accepted  views. 

The  polarization  curves  for  simple  irreversible  or  quasi-reveralble 
reactions  derived  from  the  application  of  [45]  or  [50]  will  not  differ 

3 

Other  series  expansions  such  as 

Q(p)  -  b  +  b  p  +  b  p  +  •••  +b  [62] 

0  1  2  n 

which  are  perhaps  more  obvious  than  [26]  do  not  converge  well  in  parts  of  the 
concentration  /  flux  range  and/or  generate  infinite  series  when  used  in 
Neumann's  integral.  Thus  each  term  in  the  series  [62]  itself  generates  an 
infinite  series  according  to 


I  bkJ0(aP)Pn+ldP 


^fn+2)  « 

n+i  r  I  — = —  I  >  » 

1  I  '  \  (2k+l)T(k-n/2 ) 


J 

Zk+l 


(aa) 


[63] 


The  effects  of  the  truncation  of  series  therefore  become  difficult  to  assess. 
However,  it  may  well  be  possible  to  find  series  even  more  convenient  than 
[26], 
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markedly  from  those  predicted  using  the  mass  transfer  coefficient  {18]  based 
on  the  uniform  flux  approximation^  except  In  the  high  overpotential  region 
where  the  polarization  plots  will  be  more  drawn  out  than  those  given  by  simple 
theory.  The  effect  of  including  the  potential  drop  in  solution  (using  [59]) 
will  be  to  further  draw  out  the  polarization  curves,  a  behavior  which  is  again 
in  line  with  experimental  observations:  veil  defined  limiting  currents  are 
difficult  to  observe.  It  is  likely  therefore,  that  the  major  application  of 
the  approach  outlined  here  will  be  in  the  analysis  of  data  in  the  high 
overpotential  region  and  to  the  prediction  of  polarization  curves  for  complex 
reaction  mechanisms . 

In  this  paper  we  have  applied  the  method  to  the  prediction  of  the 
behavior  of  simple  electrode  reactions  in  the  steady  state.  Numerical  data, 
applications  to  more  complex  reactions  as  well  as  tc  the  prediction  of  the 
behavior  for  a  variety  of  non-steady  state  experiments  will  be  given 


elsewhere . 
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Table 

The  coefficients  d^  (Equation  [38])  expressed  as  a  function  of  the 
coefficients  c^  (Equation  [26])  for  different  lengths  of  the  Fourier  series. 


Series 


Length 

d 

0 

d 

1 

d 

2 

d 

3 

d 

4 

n— 0 

2D(C°°-C‘)/* 

n-1 

c 

c 

0 

1 

n-2 

c  - 

c 

c 

2c 

0 

2 

1 

2 

n-3 

c  - 

c 

c  -3c 

2c 

4c 

0 

2 

1  3 

2 

3 

n— A 

c  - 

C  +C 

c  -3 

2c  -8c 

4c 

8c 

0 

2  4 

1  3 

2  4 

3 

4 

n-5 

c  - 

C  +C 

c  -3  +5c 

2c  -8c 

4c  -20c 

8c 

0 

2  4 

13  3 

2  4 

3  3 

4 

n— 6 

c  - 

C  +C 

c  -3  +5c 

2c  -8c  +18c 

4c  -20c 

8c  -48c 

0 

2  4 

13  3 

2  4  s 

3  5 

4  6 

-c 

6 

n-7 

c  - 

C  +C 

c  -3  +5c 

2c  -8c  +18c 

4c  -20c 

8c  -48c 

0 

2  4 

13  3 

2  4  6 

3  5 

4  6 

-c 

-7c 

+56c 

6 

7 

7 

n-8 

C  - 

C  +C 

c  -3  +5c 

2c  -8c  +18c 

4c  -20c 

8c  -48c 

0 

2  4 

13  5 

2  4  6 

3  3 

4  6 

-C 

6 

+C 

8 

-7c 

7 

-32c 

B 

+56c 

7 

+160c 

B 

n-9 

c  -c  +c 

0  2  4 

c  -3  +5c 

13  3 

2c  -8c  +18c 

2  4  6 

4c  -20c 

3  5 

8c  -48c 

4  6 

-c  +C 

6  B 

-7c  +9c 

7  9 

-32c 

B 

+56c  -120c 

7  9 

+160c 

8 

n-10 

C  -c  +c 

0  2  4 

c  -3  +5c 

13  5 

2c  -8c  +18c 

2  4  6 

4c  -20c 

3  3 

8c  -48c 

4  6 

-c  +c  -c 

6  8  10 

-7c  +9c 

7  9 

-32c  +50c 

8  10 

+56c  -120c 

7  9 

+160c  -400c 
8 

able  1 .  (Continued) 


16c  -112c  32c  -256c  64c  -576c  128c 

5  7  6  8  7  9  8 

+432c 


16c  -112c  32c  -256c  64c  -576c  128c  -1280c 

5  7  6  8  7  9  8 

+432c  +1120C 

9  10 


Table  2 

The  coefficients  e^  (Equation  (41])  and  the  function  f ^ (Equation  (44] 
a  function  of  j . 


Glossary  of  Symbols  Used 


I 


I 


a 

b 

l 

c 

1 

d 

1 

e 

l 

f 

l 

C 

<0 

c 

c 

Av 

c' 


D 

F 

i 

( 

1 

I 


Q 

R 

r 

s 


Radius  of  disk,  cm 

Weighting  function  series  coefficients 
Fouritf*  coefficients,  mols  (cm  s)  1 
Fourier  coefficients,  mols  (cm  s)  1 
Fourier  constant  terms 
Fourier  series  integral  terms,  cm  1 

.  3 

Concentration,  mols  cm 
Bulk  concentration,  mols  cm  3 
Average  concentration,  mols  cm  3 
Surface  concentration,  mols  cm  3 
Diffusion  coefficient,  cm2s  1 
Faraday  constant,  96485  C  equivalent*1 
Exchange  current  density,  A  cm*2 
Current,  A 

Current  density,  A  cm  2 
Bessel  functions 

Heterogeneous  rate  constant,  cm  s  1 
Heterogeneous  standard  rate  constant,  cm 

Mass  transport  coefficient,  cm  s  1 

Mass  transport  Coefficient,  cm  s  1 

Flux,  mols  cm  2s  1 
Gas  constant,  8.314  J  mols^K  1 
Radial  coordinate ,  cm 
Laplace  transform  variable 


Temperature ,  K 


Coordinate  normal  to  plane  of  disk,  cm 
Charge  of  an  ion 

Transfer  coefficient  (when  in  exponent) 

Continuous  dummy  integration  variable. 

Solution  potential,  V 
a  sin(0) ,  cm 

Solution  Conductivity,  ohms  Xcm  1 
Overpotential,  V 
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